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ABSTRACT 

Aims. Our aim is to precisely measure the physical parameters of the eclipsing binary IO Aqr and derive a distance to this sys¬ 
tem by applying a surface brightness - colour relation. Our motivation is to combine these parameters with future precise distance 
determinations from the GAIA space mission to derive precise surface brightness - colour relations for stars. 

Methods. We extensively used photometry from the Super-WASP and ASAS projects and precise radial velocities obtained from 
HARPS and CORALIE high-resolution spectra. We analysed light curves with the code JKTEBOP and radial velocity curves with 
the Wilson-Devinney program. 

Results. We found that IO Aqr is a hierarchical triple system consisting of a double-lined short-period (P = 2.37 d) spectroscopic 
binary and a low-luminosity and low-mass companion star orbiting the binary with a period of > 25000 d (> 70 yr) on a very eccentric 
orbit. We derive high-precision (better than 1%) physical parameters of the inner binary, which is composed of two slightly evolved 
main-sequence stars (F5 V-IV + F6 V-IV) with masses of M] = 1.569 ± 0.004 and M 2 = 1.655 ± 0.004 M Q and radii R t =2.19±0.02 
and R 2 = 2.49 ± 0.02 R 0 . The companion is most probably a late K-type dwarf with mass « 0.6 M 0 . The distance to the system 
resulting from applying a (V-K) surface brightness - colour relation is 255 ± 6 (stat.) ± 6 (sys.) pc, which agrees well with the 
Hipparcos value of 270 4 ^ pc, but is more precise by a factor of eight. 

Key words, binaries: spectroscopic, eclipsing - stars: fundamental parameters, distances, solar-type 


1. Introduction 


Detached eclipsing binary stars are a very important source 
of precise stellar parameters, especially of absolute radii and 
masses. It has been proven that such parameters combined with 
empirical stellar surface brightness - colour (SBC) relations 
derived from interferometr ic measurements (e.g. iKervella et alJ 
20041 IChallouf et alJl2CTl4h result in accurate distance determi¬ 
nations even to extragalactic eclipsing binaries. Indeed, very 
accurate distances to the two Magellanic Clouds determined 
with the eclipsing binary method have rec ently been reported 
(iPietrzvnski et all201 3llGraczvk et ali2014l) . In this method, the 
SBC relation is used as a tool for predicting the angular diame¬ 
ters of the system components dLacvll977l) . The distance is sim¬ 
ply derived by comparing the absolute diameter of a component 
with its predicted angular diameter. The method involves sim¬ 
ple geometric considerations and is therefore a powerful tool for 
measuring distances with a minimum of modelling assumptions. 


However, the procedure may be reversed: using prior knowl¬ 
edge of the distance to a particular detached eclipsing binary, 
SBC empirical relations can be derived in a way completely 
independent of interferometry. This basic idea has long been 


known and was firstly used by IStebbinsI dl910l 1191 ll) to derive 
the surface brightness of the two components of Algol and /3 
Aur. Within about 250 pc of the Sun there are ~ 100 confirmed 
detached eclipsing bin aries w ith paralla xes measur ed by the Hip¬ 
parcos mission (e.g. lKruszewski & Semeniukll999t (Torres et al.l 
l2010t) . However, about half of them still lack a detailed anal¬ 
ysis of their physical parameters or any analysis at all. These 
systems comprise components covering a wide range of spectral 
types from B1 to K0 that is suitable for calibrating precise SBC 
relations in the corresponding colour ranges. The are several fac¬ 
tors that influence the precision of SBC calibrations, most impor¬ 
tantly, 1) the uncertainty on absolute dimensions, especially the 
radii, 2) trigonometric parallax errors, 3) intrinsic variability of 
the star, and 4) zero-point shifts between different photometric 
systems. 

Recently, a few catalogues of precise absolute dimension 
determinations from eclipsing binar i es have been pub lished 
dTorres et alJ 1201 Ot lEker et all 1201 4t I South worth! 1201 4l) that 
cover a wide range of temperatures. However, these catalogues 
show that for the purpose of SBC calibration, a precise analysis 
of new eclipsing binaries is needed. For example, [Torres et akl 
( 2010 . see their Fig. 13) listed 15 components from nine systems 
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(El Cep, V422 Cyg, RZ Cha, GX Gem, BW Aqr, DM Vir, CD 
Tau, AI Phe, and V432 Aur) in a region of T e g - log g plane cov¬ 
ering a box of 6000 - 7000 K, 3.6 - 4.1 dex (F-type subgiants). 
Despite the number of systems with well-determined masses and 
radii in this region of parameter space, the calibration of the SBC 
for these stars is not very reliable because one system lacks a 
trigonometric parallax (V442 Cyg), three systems have very low 
quality trigonometric parallaxes (GX Gem, BW Aqr, and DM 
Vir), the Hipparcos parallax of another system (AI Phe) is incon¬ 
sistent with the value derived using the eclipsing binary method 
(lAndersen et al.!fl988l) . and one system is magnetically active 
with distortions to the light curve from star spots (V432 Aur). 
This leaves only five components in three non-active, detached 
eclipsing binaries with precise absolute dimension determina¬ 
tions within 250 pc from the Sun, and with secure trigonometric 
parallaxes. 

It is believed that parallax uncertainties from Hipparcos 
are do min ated by statistical errors - photon statistics (e.g. 
Ivan Leeuwenll2007b l. For F-type eclipsing binary with a par¬ 
allax of 4 mas we expect an error of about 1 mas, or in other 
words, a relative precision of 25 %. To improve this, we could 
average results from many components over colour bins because 
the statistical uncertainty would decrease as VfV, where N is the 
number of components used in analysis. 

The Gaia mission (iPerrvman et al.l200ll) is expected to much 
improve the Hipparcos distances by minimising statistical uncer¬ 
tainties, and for the F-type system mentioned before, the uncer¬ 
tainty of the parallax should be at most 0.5 % dde Bruiine et alJ 
l2015l l. The SBC calibration would benefit from such precise and 
accurate Gaia parallaxes obtained for a larger number of de¬ 
tached eclipsing binaries. The photometric distances derived for 
many of them would also be used as an independent check of 
GAIA results (if some significant systematic were to occur). 

Inhomogeneous photometry causes difficulties with transfor¬ 
mations between different photometric systems and zero-point 
shift uncertainties. The situation somewhat resembles problems 
occurring with the accurate determination of temperatures and 
angular diameters using the infrared flux method, where the 
main error co mes from photometric c alibration uncertainties (see 
discussion in lCasagrande et al1l2014t) . Anticipating highly accu¬ 
rate parallaxes from Gaia (precision better than 1 %) and the stan¬ 
dard precision of present-day radii determination (better then 
2%), the systematic error of future SBC calibration would be 
dominated by photometric and radiative properties of stellar at¬ 
mosphere uncertainties (Andersen, priv. communication). 

Here, we present a detailed analysis of IO Aqr (HD 196991, 
HIP 102041, a 2 ooo = 20 / '40"'45f5, <5 2 ooo = 00 o 56'2l''), an eclips¬ 
ing binary system composed of two mid-F-type stars located in 
the middle of the T e g - log g plane of the region described above. 
Despite being relatively bright (V ~ 8.8 mag), this star was not 
recognised as an eclipsing binary until it was observed with Hip¬ 
parcos and assigned its variable star name by iKazarovets et al.l 
(1 1 9991) . Time-series photometry is also available from th e Wide- 
Angle Search for Planets (WASP; iPollacco et al.ll2006h . During 
our analysis we discovered that IO Aqr is a triple system, where 
a close inner binary has a low-mass companion star on a mu ch 
wider orbit. The system was analysed before bv lDimitrov et al.l 
( 2004 ). They derived a set of fundamental parameters for the 
binary, but did not detect any trace of the tertiary component. 
Although the companion complicates the interpretation of ob¬ 
servations for IO Aqr, the companion is optically faint (less than 
1% of the total light in V) and, being on a relatively wide orbit, 
exerts only long-term perturbations on the inner system, so we 


have been able to determine precise parameters for both stars in 
the eclipsing binary component of this hierarchical triple system. 

2. Observations 

2.1. Photometry 

2.1.1. ASAS 

The V-band photometry of IO Aqr, publicly available from the 
ASAS catalogu^] (ACVS; iPoimansld! 12002) , spans from 2001 
March 27 to 2009 December 1 and contains 384 good-quality 
points (flagged “A” in the original data). This photometry was 
extended by 266 measurements from ASAS-North obtained in a 
second ASAS station on Haleakala, Hawaii. We also have access 
to previously unpublished /-band photometry from the ASAS 
and ASAS-North catalogues that spans from 1998 August 21 to 
2009 June 4 and contains 1174 good points. In total we have four 
ASAS data sets, two in V band and two in I band. 

2.1.2. WASP 

IO Aqr (1SWASP J204045.47+005621.0) is one of several mil¬ 
lion bright stars (8 ^ V £ 13) that h ave be en obse rved by 
WA SP. The WASP surv ey is described in IPollacco et al.l (120061) 
and lWilson et al.l (l2008h . The survey obtains images of the night 
sky using two arrays of eight cameras, each equipped with a 
charge-coupled device (CCD), 200-mm f/1.8 lenses and abroad- 
band filter (400-700 nm). Two observations of each target field 
are obtained every 5-10 minutes using two 30 s exposures. The 
data from this survey are automatically processed and analysed 
to identify stars with light curves that contain transit-like fea¬ 
tures that may indicate the presence of a planetary companion. 
Fight curves are generated using synthetic aperture photometry 
with an aperture radius of 48 arcsec. The data for all the stars 
in each target field are processed using the sysrem algorithm 
dTamuz et al.l 12005) to produce differential magnitudes that are 
partly corrected for systematic noise trends in the photometry. 
Data for this study were obtained between 2006 July 4 and 2010 
October 6. 

2.1.3. HIPPARCOS 

We downloaded the Hipparcos photometry in Hp magnitudes 
from a public archivcQ . We removed four discrepant high points 
from the data set and finally had 87 photometric points in total. 

2.2. Spectroscopy 
2.2.1. HARPS 

We obtained spectra of IO Aqr with the High Accuracy Radial 
velocity Planet Searcher IHARPSi lMavor et al.l2003l) on the Eu¬ 
ropean Southern Observatory 3.6 m telescope in Fa Silla, Chile. 
IO Aqr is a bright target, therefore observations were gener¬ 
ally obtained in marginal observing conditions or during twi¬ 
light. Observations were obtained between 2009 October 10 and 
2014 September 10. A total of 17 spectra were secured in high- 
efficiency (“EGGS”) mode. The exposure times were typically 
260 s, resulting in an average signal-to-noise ratio (S/N) per pixel 


1 http: //www. astrouw. edu. pi/asas/?page=aasc 

2 https://www.rssd.esa.int/index.php?project=HIPPARCOS 
&page=Epoch_Photometry 
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of ~ 60. All spectra were reduced on-site using the HARPS data 
reduction software (DRS). 

2.2.2. CORALIE 

Twenty spectra were obtained with the CORALIE spectrograph 
on the Swiss 1.2 m Euler telescope at La Silla observatory, Chile, 
between 2008 October 5 and 2010 November 1. The exposure 
times varied from 780 s to 1200 s, giving a typical S/N near 
5500A of 50 per pixel. The spectra were reduced using the auto¬ 
mated reduction pipeline. 

3. Analysis 

3.1. JKTEBOP 

We used jkteboi0 version 25 dSouthworth et al|2004allbl and ref¬ 
erences therein) to perform l east-sq uares fits of the ebop light- 
curve model dPopper & EtzelllT981l) to the various photometric 
data sets. 

For the epoch of primary eclipse for each photometric data 
set, Tq, we used the time of mid-eclipse that was closest to the 
average time of observation for that data set. Other free pa¬ 
rameters in the least-squares fit were a normalisation constant, 
the surface brightness ratio J — S i /.S' 1 , where S\ and S 2 are 
the surface brightness of the primary and the secondary star, 
respectively; the sum of the radii relative to semi-major axis, 
r\ + r2 — (R 1 + Ri)/a', the ratio of the radii, k - rilr\, the 
orbital inclination, i and the orbital period, / J (H -h- The mass ra¬ 
tio was fixed at the value derived from the spectroscopy de¬ 
scribed below. We consulted various tabulations of linear limb- 


darkening coefficients ( 

van Hamme 1993|; Dlaz-Cordoves et al. 

1995; Claret et al. 1995 

; Claret 2000) and then set the value of 


ri D for both stars to be equal to a fixed average value from these 
tabulations. We performed independent least-squares fits to in¬ 
vestigate the sensitivity of the other parameters to the assumed 
value of xld and accounted for this additional uncertainty when 
calculating the standard errors on these parameters. We assumed 
that the orbit is circular and fixed the gravity-darkening coeffi¬ 
cients of both stars to 0.2 - this parameter has a negligible effect 
on the results. 

The luminosity ratio and third light contribution estimated 
from the analysis of the spectroscopic data were included as con¬ 
straints in the least-squares fits. The relative weightings of the 
data were set equal to the root-mean-square (rms) residual of a 
preliminary fit to each data set. We used the residual permutation 
method to estimate robust standard errors on the free parameters. 
The results are given in Tabled and the fits to the individual light 
curves are shown in Fig. |T] Further details of these fits are given 
below. 

3.1.1. WASP photometry 

We divided the WASP data into eight subsets, each subset com¬ 
prising the data from one camera and one observing season. Out¬ 
liers were identified using the residuals from a preliminary light 
curve model and were removed from the analysis. The num¬ 
ber of data points used in the analysis and observation dates of 
each subset are given in Table Q] We identified between 6 and 
25 nights in each subset of data where there was an offset be¬ 
tween the data observed on that night and the rest of the data. 
We modified our version of jktebop to include an offset for the 

3 www.astro.keele.ac.uk/~jkt/codes/jktebop.html 



Fig. 1 . From bottom to top: WASP, ASAS I-band, and ASAS V-band 
photometry of IO Aqr with jktebop light-curve model fits. The numbers 
are the ranges of Julian Dates for WASP data. 


data obtained on each of these nights as additional free parame¬ 
ters in the least-squares fit. The values and standard errors given 
in Tableware the weighted mean and standard error of the mean 
for the results from these eight least-squares fits. The times of 
mid-eclipse derived from these least-squares fits are given in Ta¬ 
ble H 


3.1.2. ASAS photometry 

We fitted the four ASAS data sets independently. These data 
span a wide range of dates of observation with only one obser¬ 
vation per night. The effect of the period variations discussed 
below on the parameters derived is negligible. To derive times of 
mid-primary eclipse for the period analysis, we divided the data 
into subsets, each covering approximately 200 nights. We then 
performed least-squares fits with only Tq, J, and a scaling fac¬ 
tor as free parameters to derive the times of mid-primary eclipse 
given in Tabled 


3.1.3. Hipparcos photometry 

We used Hipparcos photometry to derive the single time of mid¬ 
primary eclipse given in Table [2] using the same method as we 
used for the ASAS data. 
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Table 3. Results of the jktebop fit to the observed LCs. Values for the WASP data are the weighted means and standard errors from 8 subsets, 
as described in Sect. 13.1.11 Symbols are defined in Sect. 13.11 Values preceded by "=" are fixed parameters in the least-squares fit. Values in 
parentheses are standard errors on the final digits of the quantity shown. Values in square brackets were constrained within the limits indicated 
during the least-squares fit. Where errors are given on fixed parameters, independent least-squares fits were performed to calculate the contribution 
of the uncertainty shown to the standard errors quoted on the other parameters. 


Parameter 

WASP 

ASAS-I 

ASAS-I-N 

ASAS-V 

ASAS-V-N 

Mean 

T 0 ( HJD-2450000) 

— 

1396.13639(24) 

3589.0060(6) 

3589.0052(6) 

5315.3727(21) 

— 

Poi-b d 

2.368145(17) 

2.3681012(51) 

2.368121(2) 

2.368123(2) 

2.368136(3) 

— 

n + r 2 

0.4238(13) 

0.4256(31) 

0.4277(56) 

0.4192(48) 

0.4243(63) 

0.4244(14) 

k = r 2 /ri 

1.1500(100) 

1.1462(39) 

1.1451(89) 

1.1551(78) 

1.1746(110) 

1.1478(36) 

i (deg) 

81.44(8) 

81.47(18) 

81.43(26) 

81.74(22) 

81.45(28) 

81.47(7) 

r\ 

0.1984(9) 

0.1938(17) 

0.1988(11) 

0.1945(25) 

0.1951(33) 

0.1976(18) 

r 2 

0.2265(18) 

0.2273(15) 

0.2268(19) 

0.2247(24) 

0.2292(33) 

0.2268(16) 

J = S 2 /S , 

0.908(6) 

0.929(5) 

0.911(4) 

0.918(12) 

0.888(18) 

— 

*LD 

=0.6 

=0.42(9) 

=0.42(9) 

=0.60(8) 

=0.60(8) 

— 

U (%) 

[1-0(5)] 

[1-0(5)] 

[1.0(5)] 

[1.0(5)] 

[1.0(5)] 

— 

latLx 

[1.21(4)] 

[1-23(3)] 

[1-23(3)] 

[1-21(3)] 

[1.21(3)] 

— 

rms (mag) 

0.014 

0.009 

0.011 

0.014 

0.018 

— 


Table 1. Log of WASP photometry for 10 Aqr. The start and end dates 
of observation, r start and r end , are JD-2450000. N ok is the number of data 
points used in this analysis, and Wiock is the number of nights of data 
that were fitted with an independent zero point in the least-squares fit. 


Table 2. Times of mid-primary eclipse for IO Aqr expressed in Helio¬ 
centric Julian Date. 
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Fig. 2. Residuals from the best-fit linear ephemeris to our measured 
times of mid-primary eclipse. 


-2485 2448161.3007+0.0010 


^start 

^end 

N ok 

^block 

-1248 

2451090.6515+0.0008 

ASAS-I 

3920 

4024 

1865 

17 

-1128 

2451374.8234+0.0002 

ASAS-I 

4622 

4753 

1339 

9 

-694 

2452402.5822+0.0013 

ASAS-I-N 

4644 

4689 

1094 

6 

-507 

2452845.4146+0.0013 

ASAS-I-N 

4970 

5117 

2251 

6 

-207 

2453555.8503+0.0014 

ASAS-I-N 

5008 

5114 

3825 

25 

-59 

2453906.3304+0.0023 

ASAS-I-N 

5008 

5114 

4944 

18 

272 

2454690.1840+0.0008 

ASAS-I-N 

5373 

5477 

5699 

22 

-779 

2452201.2893+0.0059 

ASAS-V 

5373 

5477 

4522 

25 

-657 

2452490.1977+0.0025 

ASAS-V 


-494 2452876.1983+0.0009 ASAS-V 

-193 2453589.0018+0.0016 ASAS-V 

-70 2453880.2872+ 0.0024 ASAS-V 

121 2454332.5974+0.0018 ASAS-V 

269 2454683.0784+0.0010 ASAS-V 

431 2455066.7165+0.0020 ASAS-V 

-16 2454008.1642+0.0025 ASAS-V-N 

94 2454268.6547+0.0042 ASAS-V-N 

581 2455421.9382+0.0016 ASAS-V-N 

740 2455798.4733+0.0011 ASAS-V-N 

903 2456184.4783+0.0017 ASAS-V-N 

-24 2453989.2174+0.0009 WASP 

264 2454671.2380+0.0013 WASP 

272 2454690.1845+0.0011 WASP 

422 2455045.4056+0.0007 WASP 

425 2455052.5095+ 0.0003 WASP 

427 2455057.2459+ 0.0002 WASP 

576 2455410.0988+0.0006 WASP 

577 2455412.4671+0.0005 WASP 


3.2. Period variations 

We performed a least-squares fit to data given in Table [2] to de¬ 
rive the following linear ephemeris for the time of mid-primary 
eclipse in Heliocentric Julian Date: 

HJD (UTC) = 2 454 046.0596(6) + 2.368 117 2(8) • E. 

The standard error on the final digit of each term is given in 
parentheses. The residuals from this ephemeris are shown in 


Fig. [2] It is clear that the period of this binary varies. Assuming 
that these changes are periodic, we can estimate their minimum 
period to be > 8000 d. This is discussed further in Sect. 13.61 

3.3. Radial velocities 

We used RaVeSpAn dPilecki et al .1120121) to measure the radial 
velocities of both stars in the IO Aqr binary system using the 
broadening function (BF) formalism (tRucinskilll992lll999l) . We 
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Table 4. RV measurements for 10 Aqr. Index “1” denotes the hotter 
star (primary), “2” is the cooler (secondary), and “3" is the tertiary 
component. Numbers in parenthesis give the uncertainty. BJD means 
the Barycentric Julian Date and CORAL means the CORALIE spectro¬ 
graph. 


BJD 

-2450000 

RV i 
(kms -1 ) 

rv 2 

(kms -1 ) 

rv 3 

(kms -1 ) 

Spectr. 

4745.49326 

-84.34(20) 

99.64(23) 

3.0(1.1) 

CORAL 

4746.59946 

117.50(20) 

-91.90(23) 

1.3(1.1) 

CORAL 

4746.70296 

98.85(20) 

-74.51(23) 

4.6(1.1) 

CORAL 

4747.49989 

-104.83(20) 

118.10(23) 

3.1(1.1) 

CORAL 

4747.63290 

-110.13(20) 

123.30(23) 

3.0(1.1) 

CORAL 

5086.57856 

-68.78(20) 

83.75(23) 

1.1(1.1) 

CORAL 

5086.59403 

-64.90(20) 

80.74(23) 

2.0(1.1) 

CORAL 

5087.55728 

123.67(20) 

-98.49(23) 

1.4(1.1) 

CORAL 

5087.66076 

109.02(20) 

-85.24(23) 

— 

CORAL 

5088.54850 

-107.94(20) 

120.58(23) 

1.8(1.1) 

CORAL 

5088.63650 

-110.38(20) 

122.96(23) 

2.2(1.1) 

CORAL 

5089.55659 

103.59(20) 

-79.53(23) 

2.3(1.1) 

CORAL 

5089.64296 

117.78(20) 

-93.71(23) 

0.9(1.1) 

CORAL 

5090.55337 

-37.84(20) 

54.57(23) 

— 

CORAL 

5090.62798 

-58.25(20) 

74.19(23) 

0.3(1.1) 

CORAL 

5120.57639 

129.38(15) 

-103.81(17) 

2.0(1.0) 

HARPS 

5447.66922 

99.38(15) 

-75.39(17) 

2.4(1.0) 

HARPS 

5447.75236 

79.77(15) 

-56.73(17) 

2.2(1.0) 

HARPS 

5449.50961 

103.02(15) 

-79.17(17) 

1.3(1.0) 

HARPS 

5467.49889 

-110.02(15) 

122.78(17) 

1.9(1.0) 

HARPS 

5468.50664 

112.30(15) 

-88.17(17) 

1.5(1.0) 

HARPS 

5469.51563 

-54.62(15) 

70.93(17) 

1.7(1.0) 

HARPS 

5479.51122 

-101.29(15) 

115.13(17) 

2.0(1.0) 

HARPS 

5499.54423 

128.36(20) 

-103.11(23) 

2.4(1.1) 

CORAL 

5499.57439 

127.18(20) 

-101.40(23) 

3.6(1.1) 

CORAL 

5500.53761 

-101.57(20) 

114.93(23) 

0.5(1.1) 

CORAL 

5500.57142 

-105.01(20) 

118.33(23) 

1.3(1.0) 

CORAL 

5501.53207 

85.99(20) 

-62.50(23) 

— 

CORAL 

5536.52398 

-68.88(15) 

84.30(17) 

0.7(1.0) 

HARPS 

6572.48783 

108.47(15) 

-84.14(17) 

1.6(1.0) 

HARPS 

6579.49121 

122.85(15) 

-97.65(17) 

1.9(1.0) 

HARPS 

6908.52162 

129.63(15) 

-103.69(17) 

2.4(1.0) 

HARPS 

6908.60189 

127.83(15) 

-102.24(17) 

1.2(1.0) 

HARPS 

6909.52328 

-94.03(15) 

107.85(17) 

2.1(1.0) 

HARPS 

6909.58819 

-102.62(15) 

116.23(17) 

3.0(1.0) 

HARPS 

6910.53338 

76.02(15) 

-52.86(17) 

1.4(1.0) 

HARPS 

6910.65710 

104.42(15) 

-80.37(17) 

1.2(1.0) 

HARPS 



Fig. 3. Broadening function power diagram of IO Aqr’s HARPS spec¬ 
trum taken 2010 September 29 (JD= 2455468.5) during full astronomi¬ 
cal night with the moon well below horizon. The low-temperature tem¬ 
plate (T eff = 3950 K) was used to highlight the signal from the third 
body in the spectrum - the small bump in the middle. The paraboloid 
curves are rotationally broadened fits to the BF of the primary (right) 
and the secondary (left), and the continuous vertical lines are the result¬ 
ing radial velocities of the three stars. 


used templates from the library of synthetic LTE spectra by 
ICoelho et al.l (120051) . We chose templates to match the estimated 
effective temperature and gravity of the stars in the binary. The 
abundance was assumed to be solar. We made preliminary ra¬ 
dial velocity fits separately to the HARPS and CORALIE data. 
There are very small differences in the systemic velocities of 
both components, but as they are much smaller than disper¬ 
sion of residuals, we neglected them. The resulting rms of the 
CORALIE radial velocities is about 30% larger than the rms 
of HARPS data. The line profiles of both stars are considerably 
broadened by the rotation of the stars. The projected equatorial 
rotation velocities of the primary and secondary components are 
visin; = 42.6+ 1.3 kms -1 andi’ 2 sin; = 52.3+ 1.1km -1 , re¬ 
spectively. The expected equatorial velocities for assuming syn¬ 
chronous rotation are 46.9 km s -1 and 53.6km s -1 for the pri¬ 
mary and the secondary, respectively, so it appears that rotation 
rates of both stars are very close to synchronous rotation. The 
line intensity ratio at 5500 A is 1.23. 


Table 5. Observed magnitudes of the system IO Aqr 


During this analysis, we noticed an additional faint feature 
in the BF near the radial velocity value 2 km s -1 in one of the 

Band 

total(err.) 

Ref. 

primary 

IO Aqr 
secondary 

spectra and confirmed its existence in every spectrum. The radial 
third body ve l oc it v of this feature slightlv changes with a snread of about 2 

B 

9.315(31) 

3 

10.145 

10.001 

15.92" 

km s -1 . We cross-correlated our spectra with different templates 

V 

8.834(22) 

1 

9.685 

9.506 

14.63" 

and found that the maximum strength of the feature occurs for 

R 

8.540(45) 

4,5 

9.403 

9.209 

13.81" 

an assumed template temperature of ~ 4000 K (assuming solar 

I 

8.254(19) 

1 

9.128 

8.919 

13.14* 

metallicity). The feature contributes about 1% of the total light 

J2MASS 

7.830(29) 

2 

8.726 

8.494 

12.13* 

at ~ 6000 A. Our radial velocity measurements for all three stars 

K2MASS 

7.555(33) 

2 

8.474 

8.222 

11.36* 

are given in Table |4l 


Source: 1 — this work (ASAS), 2 — ICutri et all d2003l) . 3 — iHdg et all 
d2000h (Tycho-2), 4 — iPickles & Depagnel d2010l ). 5 — IZacharias et al. I 
d2004f) (NOMAD) 
a estimated from spectra 
b assuming spectral type K6 V 


The detection of period variations (Sect. [3T2l > and third light 
in the spectra led us to the hypothesis that there is a third body in 
the system, probably a low-mass main-sequence star. This com¬ 
panion star probably is responsible for the additional signal in 
the spectrum and the period changes due to the light-time effect, 
which result from its orbital motion around the main binary. 
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3.4. Interstellar extinction 


The SIMBAD/VizieR database quotes a few different spectral 
type estimates for IO Aqr: F5 V, F5 III, F 6 V, and GO. This range 
is too broad to be useful for finding intrinsic avera ge co l ours 
of the system. Instead, we used exti nctio n maps (Sc hleg el et al. 
1998 1 with recalibration bv lSchlaflv & FinkbeineJ ( 2011 ) to de¬ 
termine the reddening in direction of IO Aqr. The total fore¬ 
ground reddening in this direction is E(B-V) = 0.070 ± 
0.001 mag. The reddening to IO Aqr is a fraction of this total. To 
estimate the reddening to IO Aqr, we assumed a simple axisym- 
metric model of an exponential disc for the distribution of dust 
within the Milky Way with scale and height lengths taken from 
iDrimmel & Spergel (1200 ll) and assumed a distance to IO Aqr of 
D = 0.25 kpc (see Sect. For Galactic coordinates of IO Aqr of 
(/, b ) = (47? 1, -23?6) and assu ming a so lar distance to the Milky 
Way centre of Rq = 8.3 kpc dGillessen et al.ll2009t) , we obtain 
E (B-VFr> Anr - 0.034 ma g . Details of our proc edure are given 
in iSuchomska et al.l (12015l) . fMaxted et ~aT1 (120141) found the scat¬ 
ter of the reddening to individual stars around the value obtained 
from these reddening maps to be of 0.034mag. 

Additionally, we used a calibration between the equivalent 
width of the int erste llar absorptio n sodium line Nal D1 and red¬ 
dening bv lMunari & Zwitteil (1 9971) . We averaged our measure¬ 
ments over several spectra of IO Aqr obtained in orbital quadra¬ 
tures between 2009 and 2014. The interstellar sodium line D1 
has one sharp component of constant radial velocity of -10.55 
km s ' 1 The mean equivalent width of this line is 0.082 A, 
which results in E(B-V)=0.019mag. The scatter of the E(B-V) 
values for individual stars around the calibration of Munari et 
al. for low reddening values is approximately 0.04 mag. Both es¬ 
timates point toward a small extinction in direction of IO Aqr. 
Finally, we assumed the extinction to be an average of both esti¬ 
mates: E(B-V)= 0.027 + 0.020 mag. 


3.5. Third light characterisation 

We estimated the colour of the third light present in the spec¬ 
trum using two spectral windows: a blue one ranging from 4000 
A to 5000 A and a red one ranging from 6000 A to 6800 A. 
These spectral windows approximately correspond to the John¬ 
son B and Johnson-Cousins R photometric bands. Each window 
was used to calculate the integrated BF profile strength of ev¬ 
ery component visible in the spectra: W \, WT, W 3 , where index 
1 is for the primary, etc. The profiles were calculated using a 
template with T e s = 6350 K, log g = 3.00, and solar metal- 
licity. We defined the third light as It, = Wil(W\ + Wi + W 3 ). 
Averaging over several spectra, we obtained Z3(B)=0.0023 and 
(3(R)=0.0078. From the observed magnitudes of IO Aqr in the B 
and R bands, we calculated the magnitudes of the third light in 
these bands (see Table0. They correspond to an observed colour 
(B-R) = 2.11 m ag and de-reddened colour (B-R)p = 2.05. 
Using Table 5 of IPecaut & Mamaiekl (|201 3lP I , we found that 
this colour corresponds to a main-sequence star of spectral type 
K 6 V. Assuming a distance of 255 pc (Sect. [4]» and a reddening 
of E(B - V) = 0.027 (Sect. 13.41 ) to IO Aqr, we calculated the 
R absolute magnitude of the companion to be M« = 6.87 mag, 
which corresponds to a spectral type between K5 V and K 6 V. 
Adding our estimate of the effective temperature given at the 
end of Sect. 13.31 we conclude that the companion has a spectral 
type K 6 ± 1 subclass. Table 0 gives the magnitude of the com¬ 


4 http://www.pas.rochester. edu/~emamaj ek/EEM_dwarf 
_UBVIJHK_colours_Teff.txt 


panion in different photom etric bands assuming colours from 
IPecaut & Mamaiekl d2013l) . The expected contribution of the 
third light in K band is ~3%. 

We found that the derived (B-R) colour depends only 
slightly on the chosen template, being on average redder for tem¬ 
plates with lower effective temperature. Assuming the same dis¬ 
tance for IO Aqr and the third light source, it is unlikely that 
the source is of spectral type later than M0 V because it would 
then be too weak to produce a significant BF signal. On the other 
hand, a star with spectral type earlier than K4 V is also unlikely 
because it would produce a much stronger BF signal that would 
match a hotter template (~ 4600 K). 


3.6. Analysis of combined light and radial velocity curves 

To perform this analysis we used v ersion 2013 of the Wilson- 
Devinney program (hereafter WD) dWilson & Deyinnevlll97ll: 
I Wilson! Il979l 19901: Ivan Hamme & WiI sonl 12007114 We cleaned 
the WASP light curve with 3-cr clipping and repeated this clean¬ 
ing three times. Then we used every third WASP observation to 
form the final light curve used in the WD analysis. This yielded 
light curves containing 4730 points in WASP band, 650 points 
in V band, 1174 points in I band, and 87 points in H p band. 

We set T, = 6475 K (see Sect. [33Tb and [Fe/H]= 0. The grid 
size was set to N - 40, and standard albedo and gravity bright¬ 
ening for convective stellar atmospheres were chosen. The stel¬ 
lar atmosphere option was used, radial velocity tidal corrections 
were automatically applied, and no flux-level-dependent weight¬ 
ing was used. We assumed a circular orbit and synchronous ro¬ 
tation for the two components. Logarithmic limb-darkening law 
were used (Klinglesmith & Sobieski 1970). The starting point 
for the parameters of the binary system were based on the so¬ 
lutions from the JKTEBOP analysis augmented by dynamical 
parameters from the preliminary solution obtained with RaVeS- 
pAn, but we also accounted for the influence of the third body 
on the orbit, as we describe below. 


3.6.1. Effective temperature estimates 


To estimate the effective temperatures of the eclipsing com¬ 
ponents, we used a number of (V - K ) and (V - I) 


colour - temperature calibrations ( 

Ramirez & Melendez! 

2001 

Gonzalez Hernandez & Bonifacio! 

2009; di Benedettol 

199a 

Casaerande et al. 20101 Masana et al. 2006: Houdashelt et al. 

2000). The intrinsic colours of the two components were calcu- 


lated from observed magnitudes, estimated reddening, flux ratios 
derived from the WD model, and from the third light contribu¬ 
tion. As a source for K-band photometry in 2MASS (ICutri et al.l 
12001 (see Table [5j. we used appropriate colour transformations 
for each calibration. Because flux ratios from the WD model de¬ 
pend on the input temperatures, we iterated the estimation of 
temperatures a few times. The resulting temperatures are aver¬ 
ages from all used calibrations. For the primary and secondary 
they are 6475 + 85 K and 6331 + 69 K, respectively, with errors 
calculated as the standard deviation. When the reddening and 
zero-point uncertainties of the photometry are added, the total 
errors of the temperature determination are 138 K and 125 K for 
the primary and the secondary, respectively. 


5 ftp://ftp.astro.ufl.edu/pub/wilson/lcdc2013/ 
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Table 6 . Results of the final WD analysis including the third body. 

Parameter 

Primary Secondary 

P (d) 

2.3681275(33) 

HJDo (d) 

2454046.0586(36) 

a(R a ) 

11.048(7) 

<7 

1.0545(14) 

i (deg) 

81.56(6) 

Q 

6.124(32) 5.681(31) 

f mean 

0.1983(13) 0.2257(15) 

y (km s -1 ) 

8.45(11) 

Teff (K) 

6475“ 6336(6) 

K (km s _1 ) 

119.80(10) 113.61(11) 

Li/Li(Hipparcos) 

1.165(15) 

L 2 /Ti(WASP) 

1.176(14) 

L 2 /Li (ASAS-V) 

1.178(14) 

Lj/LAASAS-I) 

1.213(13) 

AfHipparcos) 

0.0048" 

A (WASP) 

0.0048° 

Z 3 ( ASAS-V) 

0.0048° 

Z 3 (ASAS-I) 

0 . 0111 ° 

<?3 

0.775(76) 

m 3 (deg) 

227(16) 

«3 sin i 3 (7? 0 ) 

6354(30) 

T 03 (d) 

2449611(584) 

P 3 (d) 

30000° 

<73 =M3 sin (3 / (Mi +M 2 ) 

0.186 

RV (HARPS) rms (m s" 1 ) 

169 179 

RV (CORALIE) rms (m s* 1 ) 

209 254 

Hipparcos rms (mmag) 

14.6 

WASP rms (mmag) 

7.1 

ASAS-V rms (mmag) 

13.4 

ASAS-I rms (mmag) 

8.4 

x 2 /dof 

1.02 


" Fixed value 


3.6.2. Analysis including the third body 

We assumed that the changes i n orbital period are caused 
entirely by the light-time e ffect (IWokied Il922t llrwi a cm 
I van Hamme & Wilsonl 12007 ). Changes in the orbital period of 
the inner binary due to gravitational interaction between the 
three bodies are expected to be at least two orders of magni¬ 
tude smaller than the variation seen in Fig. l2l (Rappaport et al. I 
1201 3l . their Eq. 10). A preliminary analysis of our radial ve- 
l ocities s howed that the systemic velocity is y « 10 km s -1 . 
iDimitrov et al.l (l2004i) reported a much lower systemic velocity 
of -2 km s _1 from their observations taken in Rozen Observa¬ 
tory with the 2.0 m RCC telescope, which is equipped with a 
Coude spectrograph. Their spectra were taken in 2001, so there is 
a difference of about nine years to our spectra. We initially inter¬ 
preted the change of the systemic velocity as a sign of interaction 
between the binary and the third body. However, we quickly re¬ 
alized that the magnitude of the systemic velocity change (~ 12 
km s _1 ) is much too large: if this were real, we would expect 
much larger O-C residuals from the linear ephemeris given in 
Sect. 13.21 We concluded that there is a large systematic zero- 
point offset in velocities derived from our spectra and the Rozen 
spectra, and because their radial velocities are on average of 
lower accuracy than ours, we dropped them from our analysis. 

In the analysis we adjusted the following parameters describ¬ 
ing the inner binary of IO Aqr: the binary orbital period P, the 
zero-epoch of the primary minimum HJDo, the systemic velocity 


0 0.2 0.4 0.6 0.8 1 



Fig. 4. Radial velocity solution of IO Aqr. The lower panel shows the 
model predictions with continous lines from the WD code calculated 
for a constant orbital period of P=2.368137 d. Overplotted are radial 
velocities for the primary (denoted 1) and the secondary (denoted 2). 
The dotted line corresponds to the systemic velocity of 9.66 km s _1 
corresponding to a mean epoch of CORALIE and HARPS spectra. The 
upper panel shows residuals from the model fit. 


y, the temperature of the secondary 73, the orbital inclination i, 
the two surface potentials Qj and 0 2 , the semi-major axis a , the 
mass ratio q, and the primary star luminosity L\ in each photo¬ 
metric band. 

The O-C diagram shows that our eclipse timing and radial 
velocity data probably cover only a small fraction of the orbital 
period of the third body, P 3 . In this case, the third body orbital 
parameters are not uniquely determined. A striking feature of our 
radial velocities is the lack of significant changes in the systemic 
velocity of the inner binary or the mean radial velocity of the 
third body. Moreover, the shape of the O-C diagram suggests 
a highly eccentric and long-period third body orbit. It was also 
clear from our trial solutions that the orbital period P 3 and the 
eccentricity e 3 strongly correlate, as is the case for the longitude 
of periastron w 3 and the time of the upper spectroscopic con¬ 
junction T03. For this reason, we decided to find solutions within 
a grid of possible periods P 3 . The additional fitted parameters 
were those describing the third body orbit: eccentricity e 3, lon¬ 
gitude of the periastron M 3 , semi-major axis <23 , and the moment 
of the upper spectroscopic conjunction T03. We also assumed 
that the outer orbit is seen edge-on (z '3 = 90 deg). 

We decided to divide the third body parameters into three 
groups - {e 3 , T 03 }, { 013 }, and {a^, T 03 }. We iteratively fitted 
these three parameter groups with the parameters of the inner 
binary (P, HIDq, y, 73, i, Q|, Q 2 , a, q, and L\ ) until we achieved 
no improvement in y 2 for a given value of P]. We repeated this 
procedure for many orbital periods P 3 within a range of from 
8000 to 50000 days. We could not find any acceptable third body 
solutions for periods shorter than about 25000 days. The reason 
is that the radial velocity difference between the inner binary 
systemic velocity and the radial velocity of the third body (AV r ~ 
8 km s -1 ) causes a tension with the shape of the O-C curve. 
The tension is only relaxed for P 3 > 25000 days. Periods longer 
than 50000 d require a very high eccentricity of the outer orbit 
(e 3 > 0.85). In all cases the inferred minimum mass of the third 
body is about 0.6 M 0 . The specific solution for P 3 = 30000 
days is reported in Table[ 6 ] Figure[4]presents the fit to the radial 
velocity data given in Table[4]assuming a constant orbital period 
during the time span of our spectroscopic data. 
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Fig. 5. Left panel: O-C residuals of the observed times of primary minimum of IO Aqr assuming the linear ephemeris from Tableland predicted 
O-C caused by the light-time effect from our third body solution - continous line. Right panel: predicted radial velocity changes of the binary 
mass centre (continuous line denoted 1+2) and the third body (dashed line denoted 3). The horizontal line corresponds to the heliocentric frame 
systemic velocity of the whole triple system. The individual binary mass centre velocities from HARPS are marked by filled circles and CORALIE 
velocities by open circles. The velocities of the third body were binned into seasonal means. 


Figure [5] presents the fit to the O-C residuals of observed 
minima times from Table [2] and the expected radial velocity 
changes of the inner binary mass centre of IO Aqr and the third 
body. To calculate the binary centre-of-mass velocities Vi +2 , we 
used radial velocities from Tableland the equation 


Vi+2 


Vi + qVi 

1 + q 


( 1 ) 


where V\p denotes radial velocities for the primary and the sec¬ 
ondary, respectively, and q is the mass ratio. The asymmetric 
minimum in the O-C diagram coincides with the periastron pas¬ 
sage of the third body near JD 2452600. The O-C diagram 
shows that our final solution for the third body orbit does not 
fit the Hipparcos data. Forcing a good fit for this one outlier re¬ 
sults in stronger binary centre-of-mass radial velocity changes 
and produces significant residua in the radial velocities. A possi¬ 
ble explanation is that we underestimated the error on Hipparcos 
minimum timing. It is also clear that we only cover a relatively 
small part of the whole third body orbit, and much room for im¬ 
provement probably still exists. 



Fig. 6. H-R diagram with the primary (1) and the secondary (2) com¬ 
ponents of IO Aqr marked. Overplotted is the PARSEC isochrone that 
provides the best match to M, L, and r cff of the two stars, corresponding 
to an age of 1.88 Gyr. Dots on the isochrone correspond to the position 
of the best fit. The zero-age main-sequence (ZAMS) is also shown. 


3. 7 . Physical parameters 

The physical parameters of the components of IO Aqr are 
presented in Table [7] Our spectral type estimate is based 
on the derived effective temperatures and the calibration by 
iPecaut & Mamaiekl(l201 3i): it is co n siste nt with the spectral clas¬ 
sification given by iHouk & Swift) d!999t) . The mean radii given 
are the radius of a sphere with the same volume as the star. The 
distortions from a perfect sphere defined as (r po i nt - r po i e )/r mean 
are 2.3% and 3.2% for the primary and the secondary star, re¬ 
spectively. Errors on luminosities were calculated from uncer¬ 
tainties of the distance estimate and the reddening and bolomet- 
ric corrections used. T he bolometric corrections reported in Ta¬ 
ble [7] were taken from lCasagrande et al.l (120101) . The fractional 
precision of derived masses is 0.25% and 0.24% for the primary 
and th e seco ndary, respec tively. 

[Dimitrov et al.l J2004 ) reported two possible solutions for ab¬ 
solute parameters of IO Aqr denoted as A and B. It becomes ev¬ 
ident that our solution lies somewhere between these solutions. 
However, we were able to significantly refine the parameters of 


the system. The radii of the two stars are too large for their mass 
or for their temperature if we were to assume that they are ma in- 
sequence stars, as was also pointed out bv lDimitrov et al.ld2004l) . 
We compared our results with stellar isochrones from th e Pa dova 
and Trieste Stellar Evolutionary Code (PARSEC; iBressan et al.l 
120121 FI We fitted masses, luminosities, and temperatures of the 
components. As we have no confident metallicity estimate for 
this star (the SIMBAD database quotes some estimates made un¬ 
der the assumption that IO Aqr is a single star), we also fitted this 
parameter. The position of the stars in the Hertzsprung-Russell 
(H-R) diagram and their measured masses are only matched us¬ 
ing isochrones for metallicities [M/H] ~ +0.3 - see Fig. [6] The 
estimated age of the sys t em is 1 .88 + 0.10 Gyr. We confirm the 
finding of lDimitrov et al.l (2004) that the two stars are still before 
the turn-off point in their evolution. 

Our absolute parameters do not depend on the orbital so¬ 
lution for the third body. The third body only affects the tim- 


6 http://stev.oapd.inaf.it/cgi-bin/crad 
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Table 7. Physical parameters of IO Aqr components 


Parameter 

Primary 

Secondary 

Tertiary 

Spectral type 

F5 V-IV 

F6 V-IV 

~ K6 V 

M (M q ) 

1.569(4) 

1.655(4) 

0.6(1) 

R(R e ) 

2.191(15) 

2.493(17) 

— 

logg (cgs) 

3.952(6) 

3.863(6) 

— 

Teff (K) 

6475(138) 

6336(125) 

-4000 

L (T 0 ) 

7.59(55) 

9.00(65) 

— 

M v (mag) 

2.57(8) 

2.39(8) 

- 7.7 

BCy (mag) 

-0.021 

-0.030 

— 

Distance (pc) 


255(9) 


Distance modulus (mag) 


7.030(72) 


Parallax (mas) 


3.93(13) 


E(B-V) (mag) 


0.027(20) 



ing of the eclipses and, to a lesser extent, the estimates of the 
distance and effective temperatures. Our solutions without the 
third body (Sect. 13.1.1b yield model parameters that are consis¬ 
tent with those reported in Table [6] The high precision of the 
derived physical parameters makes IO Aqr another eclipsing bi¬ 
nary star on a list of systems that have a precision of the abso lute 
dimension determination better than 3% (e.g. iTorres et al.l2010h . 


4. Distance and space velocity 

The distance to the system was derived using two calibra¬ 
tions of visual surface brightness versus (V — K) colour rela¬ 
tions dKervella et al.l 120041 Idi Benedettoll2005i) based on mea¬ 
surements of the interferometric stellar angular diameters. The 
2MASS K -band magnitude was tranformed onto the John- 
son photometri c syste m u sing transformati on equations from 
iBessell & Brett! (119881) and ICarnenterl (|2001|PI The magnitudes 
were de-reddened with E(B — V) = 0.027 mag. To derive indi¬ 
vidual K-band magnitudes of the components, we used the K- 
band light ratio extrapolated from the WD model h\(K) = 1.262 
and the third light contribution h{K) - 0.030. The resulting de- 
reddened Johnson V-band and K-band magnitudes for the two 
stars are Vi = 9.601 mag, Vi = 9.422 mag, and Ki = 8.505 mag, 
K 2 = 8.252 mag. The distance to IO Aqr der ived f rom the SBC 
relatio ns reported by iKervella et al.l (120041) and Idi Benedetto! 
d2005l) are 257.2 pc and 252.2 pc, respectively. The resulting av¬ 
erage distance is 255 pc, corresponding to a distance modulus 
m — M = 7.030 mag. 

The systematic errors contribute as follows: the uncertainty 
of the empirical calibrations of the surface brightness (0.040 mag 
error in distance modulus), the uncertainty of the extinction law 
(0.005 mag), the metallicity dependence of the surface bright¬ 
ness relations (0.004 mag), and the zero-point errors of V- and 
K-band photometry (0.030 mag in total). When we combine 
the contributions in quadrature, we have a total systematic er¬ 
ror of 0.051 mag. The total statistical error comes from the un¬ 
certainty in the semi-major axis (0.005 mag), the uncertainty 
of the sum of fractional radii r\ + n (0.010 mag), the uncer¬ 
tainty of the third light (0.022 mag), the errors of V- and K- 
band magnitudes (0.007 and 0.033 mag, respectively), the red¬ 
dening uncertainty (0.015 mag), the error of the mean of the two 
calibrations (0.021 mag), and from separating the magnitudes 
(0.002 mag). When we combine these uncertainties in quadra¬ 
ture, we have a total statistical error of 0.049 mag. When the er¬ 


rors in distance modulus are translated into errors in parsecs, the 
final result is 255 ± 6 (stat.) ± 6 (syst.) pc. The distance corre¬ 
sponds to a parallax of 3.93 ± 0.13 mas (total error). 

The important consistency check of model parameters we 
derived (especially their radii and luminosity ratio in different 
bands) is the distance to each of the system components re¬ 
turned by our procedure. In our case, the distance moduli of the 
two stars agree excellently well: they are different by less than 
0.001 mag for each of the surface brightness calibration used. 

The parallax to the binary was determined by the Hippar- 
cos sat ellite mission. The first published full data reduction 
dPerrvman et al.l 1 19971) gave 5.42 + 1.46 mas. Subsequently, a 
new data reduction ( van Leeuwenll2007allbh reported a dist inctly 
smaller parallax of 3.71 ± 0.94mas. iDimitrov et al.l (120041) used 
an eclipsing binary bolometric flux scaling method and derived 
a distance to IO Aqr of 256 pc, corresponding to a parallax of 
3.91 mas. However, they neither included any third light in their 
considerations nor gave a distance error estimate. Both of the 
latter distance determinations agree very well within the uncer¬ 
tainties with our derived distance._ 

The proper motion of IO Aqr given bv lvan Leeuwenl d2007bl) 
is (yU [V cos 6,/-is) = (-4.59 ± 1.00,-22.42 ± 0.84) mas yr _1 . 
That proper motion expressed in Galactic coordinates is 
(/// cos b, /Jb) = (— 21. 7 ± 1.2, -7.4±0.4)mas yr~* using the trans¬ 
formations given by iPoleskii 1 20131) . The transverse velocity in 
Galactic coordinates is (-26.2 ± 1.5,-8.9 ± 0.5) km s -1 . He¬ 
liocentric Galactic space velocity components, that is, not cor¬ 
rected for solar peculi ar moti on, were calculated using the equa¬ 
tions given by iJohnson & Soderbloml (119871) . and we obtained 
(u,v,w)= (21.8 ± 0.6, -14.9 ± 1.6, -11.5 ± 0.8) km s _1 . The val¬ 
ues and uncertainties given do not take into account possible sys¬ 
tematic errors caused by a long-term proper motion drift of the 
inner binary induced by the third body. At the moment, we have 
insufficient information to correct for this effect. 


5. Final remarks 

With the high-precision data now available, we have been able to 
show that IO Aqr is a hierarchical triple star system with a low- 
luminosity and low-mass companion on a wide eccentric orbit. 
For the purpose of accurate SBC calibration, pure eclipsing bi¬ 
naries (i.e. without stellar companions) are preferred. However, 
when the physical properties, especially the luminosity, of a ter¬ 
tiary companion are well determined or/and the flux of the com¬ 
panion in optical and NIR is negligible, a given eclipsing binary 
can be used safely. A criterion is the degree of calibration preci¬ 
sion we wish to obtain: 1% of precision demands that the tertiary 
contribution in V and K bands should be at most comparable 
with this limit. The companion of IO Aqr is faint enough in the 
optical to be ignored, and its estimated K-band contribution is 
still small, so that it has a negligible effect on the calculated an¬ 
gular diameters of the IO Aqr components. Using the SBC cal¬ 
ibration of ldi Benedettol (l2005h . the predicted angular diameters 
are 0.081 ± 0.002mas and 0.092 ± 0.002 mas for the primary and 
the secondary, respectively. 

We derived high-precision absolute dimensions of IO Aqr 
that are among the most precise parameters of F-type stars pub¬ 
lished to date. Nonetheless, our temperature determination (pre¬ 
cision 2%) can be improved by using dedicated Stromgren pho¬ 
tometry or by a detailed atmospheric analysis of the separated 
spectra. We leave this to future work on IO Aqr. 


7 http://www.astro.caltech.edu/~jmc/2mass/v3/ 
transformations/ 
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